Photon correlations in positron annihilation 
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The two-photon positron annihilation density matrix is found to separate into a diagonal center 
of energy factor implying maximally entangled momenta, and a relative factor describing decay. 
For unknown positron injection time, the distribution of the difference in photon arrival times is a 
double exponential at the para-Ps decay rate, consistent with experiment (V. D. Irby, Meas. Sci. 
Technol. 15, 1799 (2004)). 
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I. INTRODUCTION 

When an electron and a positron with opposite spin 
annihilate, two correlated photons with total energy 2 x 
.511 MeV are created. These annihilation 7-rays cannot 
be manipulated using optical beam splitters and mirrors, 
so interference experiments and applications in quantum 
information are not practical. However, positron annihi- 
lation is important in medicine and material science 1]. 
In medical imaging, coincident detection of the annihi- 
lation photons is the basis for positron emission tomog- 
raphy (PET). In material science positron annihilation 
spectroscopy (PAS) gives information on electron den- 
sity and the distribution of electron momenta. 

Positrons are created by the decay of radioactive nuclei 
such as '^'^Na or ^^F imbedded in the sample of inter- 
est. For example, the 1.275 MeV nuclear 7-ray emitted 
immediately following the positron emission from '^'^Na 
determines the time of positron injection. In positron 
lifetime (PAL) measurements the arrival time difference 
between the nuclear photon and one of the annihilation 
photons is measured. Positron annihilation in condensed 
matter proceeds through bound states of positrons with 
electrons, atoms, molecules and various defects [l[. The 
annihilating positron and electron form a free or bound 
hydrogen-like positronium (Ps) atom. In vacuum, sin- 
glet or para-Ps decays into two 7-rays with a lifetime of 
125 ps. In a— Si02 the para-Ps lifetime is increased to 
156 ps due to modification of the dielectric constant and 
electron mass relative to vacuum 0. 

Recently it has been suggested that measurement of 
the arrival time difference between paired annihilation 
photons will improve signal to noise in medical imag- 
ing applications, leading to time of flight (TOF) PET 
3. This is plausible because the most widely accepted 
viewpoint is that the minimum quantum uncertainty in 
time is zero due to detection-induced nonlocal collapse 
Q. Irby measured the time interval between detection 
of the annihilation photons from a ^^TVa source and ob- 
tained 123 ± 22 ps [1]. This is a surprising result since, 
in his experiment, the annihilation photons originate in a 
source a few mm thick and a photon travels almost 4 cm 
in air in this time. 



'Electronic address: ' margaret . hawton@lakeheadu. ca] 



To explain these observations, Irby generalized the 
Einstein, Podolsky and Rosen (EPR) Q example of po- 
sition and momentum as elements of reality to include 
time and energy dependence @ . Using entangled spins as 
an illustration, he showed that restriction of one observ- 
able leads to reduced nonlocality of its conjugate. He at- 
tributed his experimental results to maximally restricted 
photon momenta, leading to elimination of nonlocality 
in the conjugate position observables. However, a com- 
plete explanation requires a theory of the 123 ps wide 
distribution of time differences that he observed. Here 
we give a quantitative explanation of his observations by 
performing a detailed analysis of Ps decay. 



II. THEORY 

This section is based on Sakurai's theory of positron 
annihilation (t}, summarized in Subsection A, trans- 
formed to relative and center of energy coordinates in 
Subsection B, and modified to explicitly include expo- 
nential decay in Subsection C. Natural units in which 
h — c — 1 are used, the electron/positron mass is de- 
noted as m, and the positron charge is e. The dimension- 
less fine structure constant is then a = e^/47r — 1/137. 
The subscript + refers to the positron and — to an elec- 
tron. We consider a relativistic expansion in powers of 
the Fermion speeds, (3+ and /3_, denoted /3± where, to 
first order in j3±, the annihilation photons are counter- 
propagating. To simplify the equations it is assumed that 
the photon pulses are well separated from the positron 
source when they reach the detectors. 



A. Positron annihilation 

Position annihilation according to the Dirac equation is 
discussed by Sakurai. He performs a perturbation expan- 
sion in powers of e and finds that the first nonzero term is 
of second order. The Feynman diagram of such a process 
is sketched in Fig. 1: An electron with four- momentum 
p_ — {E-,p-) is scattered to four-momentum q — 
((7o,q) at space-time point X2 — (ti,yiri) while emitting 
a photon with four- momentum fc2 — (a;2,k2) . At x\ this 
electron annihilates with the positron and emits a photon 
with four-momentum k\. If instead the positron is scat- 
tered first, q —q and the photons are interchanged. 
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B. Relative and center coordinates 



ki 
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FIG. 1: A two photon Feynman diagram. An electron, p-, 
emits a piioton, k2, wliile scattering to a virtual state, q. It 
then annihilates with a positron, p+, while creating a second 
photon, ki. 



Sakurai obtains a scattering cross section for two-photon 
annihilation of 7rrQ//3+ where tq = a/m. The Bohr ra- 
dius, uq = 1/ (am), is larger than rg by a factor a~^, 
so the volume of an atom appears infinite on the length 
scale rg and the center of energy momentum is conserved, 
that is 



ki + k2 = P4 



(1) 



Sakurai applies his scattering theory to Ps by setting 



the electron density equal to {ipis 
obtains a decay rate 



1/ 7r(2ao)' 



r 



1 



and 



(2) 



equivalent to a lifetime F ^ = 125 ps. 

In Sakurai's covariant formulation energy and momen- 
tum are conserved at the vertices and the state q de- 
scribes a virtual particle for which the Fermion disper- 
sion relation is not imposed. However, since the final 
and initial states describe real particles, the dispersion 
relations 



= mi , 



(3) 



E+ = A/m2 + |p±|' 



must be satisfied. In the more usual noncovariant for- 
mulation of perturbation theory, the dispersion relation 
is satisfied by the virtual Fermion but energy is not con- 
served between ti and t2. To zero order in f3± the an- 
nihilation photon fc2 has energy to, so the excess en- 
ergy of the virtual state must be greater than m. Thus 
the intermediate state in Fig. 1 persists for less than 
TO~^ = 1.3 X 10~2^ s, implying that two photon annihi- 
lation is effectively instantaneous. 



Here the center (of energy) and relative coordinates. 



kc = ki+ k2, = i (ki - k2) , 



(4) 



P , Pr 



(P+ -P-), 



Xc = - (xi -I- X2) , and x,. = xi - X2, 

will be used. Since ki • Xi + k2 • X2 = kc • Xc + kr • x^ 
for the photons and p+ • xi + p_ • X2 = Pc • Xc + Pr • x^ 
for the Fermions, the exponent in a Fourier transform 
is preserved by this transformation, and relative momen- 
tum and position are conjugate observables, as are center 
momentum and position. 

For counterpropagating photons the magnitudes of ki 
and k2 should be added (subtracted) to obtain the mag- 
nitude of the relative (center) wave vector so that, ac- 
cording to ^ and dU, 



UJ : 



0J2 



2 Ik, 



To second order in f3± the Ps total energy is 
E = 2m + pi/Am. 



(5) 



(6) 



For a positron created at time Iq contributions with 
different pc rapidly get out of phase due to the factor 
exp [—iPc {t — to) /4to] , leading to a density matrix that 
is diagonal in center of energy momentum. The relative 
dynamics, described by k,, is decoupled from the center 
motion, described by k^. In relative and center coordi- 
nates conservation of momentum, ([1]), becomes 



kc = Pc 



(7) 



Since Pc has a definite value, the momenta of the an- 
nihilation photons are maximally restricted according 
k2 Pc — ki as observed by Irby. 



C. Dynamics 

Sakurai calculates the Ps decay rate so, implicitly, uj 
isn't exactly equal to E, but has a linewidth, F. Decay 
as a function of t will be considered in this subsection. 

A pure state will be written as a linear combination of 
a Ps atom in the Is state with definite center of mass mo- 
mentum Pc, and the two annihilation photons described 
by their relative and center momenta. If a positron is 
injected at time to the Schrodinger picture (SP) state 
vector is then 



|*kj =Cis (t)|ls,ke) + 



(i)|k„ke). (8) 



for t > to and l^'k^) = for t < io- We will take the 
volume, V, to be finite so that the momenta are discrete. 
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To second order in e the dynamical equations describing 
the relative motion for t > are ^ 



Ck. (t) = -iwck,, (t) - zC/,(2) ci, (i) 



(9) 



where the dot denotes differentiation with respect to t 

■ (2) 

and U = Ur S (kc — Pc) is the time derivative of 
the transition matrix element from Ps to the two-photon 
state. Eqs. (O describe Weisskopf-Wigner spontaneous 
emission that is exponential in time and Lorentzian in 
frequency. A system of equations of the form © are 
solved in the interaction picture in For t — ^ 
decay is essentially complete so that the photon pulse is 
well separated from the source and [l0[ gives 



Ck 



exp [— icj (t — io)] 
Lj- E + iT 



(10) 



in the SP with w = 2 |kj.| and E — 2m to first order in /3±. 
The factor A is a constant and ([TU]) can be normalized 
using the integral II in Appendix A with the result 



SttP cxp [— zo; {t — to)] 



VE^ Lo- E + iT ■ 
A pure state vector is of the form 

|*kj = e (ri - to) e (t2 - to) |k,) (g, 



(11) 



(12) 



Substitution of k = k^ , r = and t = t ~ to in integral 
12 of in Appendix B gives 



V'(|x.|,i) = 




(17) 



X exp 



{iE + r)[t-to' - |x. 



where a similar term involving t — to + ^\xr\ has been 
neglected. This wave function is normalized if it is as- 
sumed that the photon pulse has propagated far enough 
so that exp [-P {t - to)] < 1. 

For a measurement described by the operator O, the 
expected value is 



kc 



o 



(18) 



where jvE'ke) given by is a pure state and the proba- 
bility for center of mass momentum kc is . Normaliza- 
tion is such that (xj.|x^) = {xr — xj,) , (kclk^) = (^kc,k'^ 
and Lke^kc = LkJckJ^ = J d^Xr\ip {xr,t)f = 1. The 
0-functions in limit the volume that the j*^ photon 
can occupy to V ~ (t — to)^ . For finite volume con- 
servation of momentum, ([7]), is approximate, with uncer- 
tainty of order tt / (t — to) in each of its components. 



III. APPLICATION TO EXPERIMENTS 



t - IXi 



(13) 



where Xj is the position of the j photon, Tj is its emis- 
sion time, the 0-functions ensure that no photons exist 
before the positron is injected, and 



In this Section, Eq. (flSj) will be applied to Doppler 
broadening (PAS experiments) and the arrival time dif- 
ference between the nuclear photon and one of the anni- 
hilation photons (PAL experiments) , and the Irby exper- 
iment will be analyzed. 



SttP ■r-^ exp [—iuj it — to)] 



VE^ 



u~E + iT 



describes the relative dynamics. 

The space-time wave function is ip{'x.r,t) 
such that 



with 



ip (Xr,t) 



cfixrip (Xr, t) |x,,) 



' 4r 6xp (iujto) 
E^^ Lj- E + iT 

k. 



exp (ik,- • Xr — iujt) 
372 ■ 



(14) 

(Xr|*r) 
(15) 

(16) 



A. Doppler broadening 

Ref. [m reports measurement of the distribution of 
the Ps center of mass momentum, so that O = |kc) (kc|. 
Substitution in (|18l) gives the probability of center wave 
vector kr as 



(19) 



This experiment was performed using a positron source 
embedded in biological tissue, and the Gaussian distri- 
bution 



P(kc 



1 



(2^) 

Strictly speaking, the k,.-amplitudes should be weighted 

as in a Is state, but F ^ Aq"^, so this can be ignored. P (k) — pi^^V/ {2Tr)" 



(7ra2) 



2^3/2 



cxp 



(20) 



with (7 = 2AkeV = 0.005m was obtained. The continu- 
ous distribution in related to the discrete probability by 
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If these center of mass momenta were to add coher- 
ently, the time uncertainty for the second photodetec- 
tion event would be very small. However, the photon 
momenta are maximally correlated so, if Xc were to be 
measured, (fT8)) gives 



(|x,) (x,|)=^pkj(x.|k,)|^^ 

k, 



(21) 



This implies that the photon center of energy is equally 
likely to be found anywhere within the allowed volume, 
since the only information available about its position is 
a consequence of causality and knowledge of the position 
and time of positron injection. 



B. PAL experiments 



C. Irby experiment 

In the Irby experiment, illustrated in Fig. 2, photons 
are emitted by a source, S, approximately imm thick. 
They are detected at the fixed positions Xi and X2 as a 
function of ii — t2 where tj is the time when a photon is 
counted at detector j. Irby derived a wave function that 
generalizes the example considered by Einstein, Podol- 
sky and Rosen (EPR) by including time dependence and 
conservation of energy 6]. He assumed zero center of 
mass motion so that the photons have momentum p and 
—p. The relative position, x^, corresponds to xi —X2 and 
the Fourier amplitude, Ck^, given by (jlip corresponds to 
f{p) in Irby's Eq. (13). 

Following EPR and Irby [5, 6] and using in the 
form |xr| = |x<| + |x>|, the wave function (fT7|) can be 
written as 



In PAL experiments such as the measurement of 
positron lifetime in a-Si02 photons are counted at 
fixed xi as a function t — tQ. It is assumed here that para- 
Ps forms as soon as the positron is injected, although in 
reality the situation is more complicated than this. To 
first order in /3± the wave vector has length m and 
arbitrary direction. The wave vector kc has a definite 
value and its magnitude is distributed according to ([20l) . 
Substitution of 6 = |xi) (xi| , T = / (fx2 |x2) (xal, ([T2|) . 
and (HZl) in (HH) gives 



(|xi)(xi|) 



r 



exp [-r {t - to)] 



X 8 (i — to — |xi 



(22) 



X J d^X2e{t-to-\x2\) 

X Ix^r^exp [-P (i - to - |xr|)] . 

This is just the trace of the density matrix over the unob- 
served second photon. If the z-axis is chosen parallel to 
ki , the distribution of k2 values is centered aX 9 = tt and 
the factor exp (P jx^l) selects solid angle determined by 
r and centered about cos 6* = — 1. To first order in /3± 



X2 



(23) 



In the limit |xi| 3> P^^, consistent with our assumption 
that the pulse is well separated from the source, jx^l « 
2 |x2| and the probability density to count a photon at 
xi a time t — to after positron injection reduces to 



(|xi)(xi|) 



exp[-P(t-to-|xi|)] (24) 



IQttV 

X e(t-to- |xi|) 



V'(|x,|,t)= / da;(5(|x<|~a;)V(|x>|+2;,t) (25) 
Jo 

where (5(|x<| — x) is a positon eigenvector with eigen- 
value X, x< is the position while t< is the time of the 
first photodetection event, and x> is the position of the 
second photon. When the first photon is counted at time 
t< the wave function collapses to the coefficient of the 



(5-function in (|25p . To ensure propagation at the speed 
of light this one-photon exponentially decaying pulse can 
be written as 



^(|x>|,t) 




(26) 



|x<|) 



--{lE + T) (t-to-|x>|) 



Time and distance dependence for the undetected photon 
is described by the last exponential, so the probability 
density is proportional to exp [— P (t — to — |x> |)] or zero. 
If the second photon is counted at time t> , allowing for 
the Xc density the probability density for coincident 
photodetection is 



1 

P = — 
V 



^|x,.| , — 



tl+t2 



(27) 



where |xr| is the detector separation, t< -f t> = ti + t2 
and ip is given by ([T71) . 

Essentially the same result is obtained from the second 
order Glauber correlation function [l2i. 



where V — {t — to)^ ■ Thus the rate at which corre- 
lated nuclear and annihilation photons are counted de- 
cays exponentially. The coefficient of the exponential 
refiects our limited knowledge of the position of the two- 
photon center of energy. 



G(2) (xi,X2) = (S(-) {xi)E^-^ {x2)E^ + '> {X2)E^ + ^ (xi)) 

(28) 



where 



(tj,Xj). For photodetection at times ti and 



t2, the positive frequency electric field operators in G'-^-' 
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d: 




^1 




FIG. 2: Irby experiment. A positron is created in the source, 
S, and the time difference between annihilation photons ar- 
riving at detectors dl and d2 is measured. 




result in a factor 

exp [—i {ujiti + ^2^2)] = cxp 



tl+t2 



FIG. 3: Comparison of exponential of the absolute value, 
(2r)~^exp(-r|a;|) with a Lorentzian, (ttF^)"^ {x^ +r^)~^ , 
(29) and a Gaussian, (Ftt)"^/^ exp (-r^a;^). 



Aw , 
+ — (ii-i2) 



Since ^yuJiuj2 = m is a constant to first order in /3± , 

2 



(30) 



equal to P given by ([27]). 

The probability density P is proportional to 
exp [— r {ti +t2 — 2io)], but Irby measured the distribu- 
tion of ti — t2, and neithe r (|27p nor the absolute square 
or Irby's wave function in [J| gives their probabilities di- 
rectly. The resolution to this problem lies in averaging 
over the positron injection time, to, that is not measured 
but must be earlier than both ti and T2. If it is assumed 
that positrons are injected at a constant rate r = 1/T , 
substitution of pT|) in (|?7)) gives 



rT 



T/2 



dtn 



(31) 



2 / L^^O 
47r|Xr| V 

X exp [-r (h +t2- 2to - |x^|)] 
X e (ti - to) e {t2 - to) . 

The integral (PT|) is evaluated as /3 in Appendix C with 
the upper limit of the to integral is taken to be the earlier 
photon emission time. The result is 



P 



exp (-r In -T2I), 



(32) 



where Tj = tj — |xj|. 

Irby fit his date to a Lorentzian curve while, accord- 
ing to ([5^ . the experimental picosecond timing analyzer 
(PTA) spectrum in Figs. 4 and 5 of Ref. '3| is a dou- 
ble exponential. This discrepancy is addressed in Fig. 
3 that shows a comparison of a double exponential to a 
Lorentzian and a Gaussian. The double exponential gives 
the sharp peaks observed by Irby while behaving like the 
Lorentzian that he used in his fits in the tails. The Gaus- 
sian has an appreciably different shape and does not fit 
the data as noted by Irby. Eq. ([5^ derived here should 
give an improved description of the experimental results. 



IV. CONCLUSION 

This paragraph describes the details of the present cal- 
culation in relation to the previous theoretical work: In 
Refs. Q and [6| the center of energy momentum is set 
equal to zero, the wave function is given as a function 
of the relative coordinates, and the time during which 
the photons interact, here to to to + F^^, is assumed to 
be known. In the present calculation the momentum of 
the center of energy has a wide range of definite values 
consistent with the PAS experiments, and the positron 
injection time is unknown. In Q all relative momenta 
are given equal weight. Since the time when the parti- 
cles interact is known, when one of the counterpropagat- 
ing photons is detected the position of the second photon 
is determined exactly and nonlocally by collapse of the 
wave function. Here and in Q the relative momenta, 
p — |kr|, are restricted by a function / (p) which we find 
here is a Lorentzian with center at |kj.| = m and FWHM 
2F, resulting in exponential decay in space-time. 

Irby attributed the unexpectedly wide range of anni- 
hilation photon PTA detection time differences that he 
observed to maximally restricted photon momenta, lead- 
ing to the elimination of nonlocality in the conjugate po- 
sition observables [6j. Here the pure states have definite 
center of energy momentum and Ps decay is described in 
terms of the relative coordinates. After averaging over 
the rmobserved positron injection time, the annihilation 
photon coincidence rate was found to be proportional to 
exp (— F |ti — T2I) where Tj is the photon emission time. 
This supports Irby's observation [4| that annihilation 
photon pulse width is limited by the Ps lifetime. Only 
the peak of the double exponential function is determined 
by the position of the positron source. This is counter to 
expectations, and should be taken into account in TOF 
PET imaging. 

Annihilation photons have played a significant role in 
the development of our understanding of quantum corre- 
lations. Their polarization correlations were considered, 
and discarded, as a candidate for the first experimentally 
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realizable test of Bell's theorem [l3|. EPR used posi- 
tion correlations of a pair of counterpropagating particles 
as their primary example of nonlocal collapse [5|. Irby 
performed a direct measurement of annihilation photon 
space-time correlations and concluded that their nonlo- 
cality is erased by maximal restriction of their momenta. 
Here we find that their momenta are maximally corre- 
lated because their center of energy momentum has a 
well defined value. Position entanglement is ascribed to 
the relative coordinates, augmented by causality. The 
observed 123 ps pulse width is attributed to uncertainty 
in the time of photon pair creation due to Ps annihilation. 



Appendix C: Irby experiment fo-integral 



We need 



rT/2 

13= dto exp [-r (2<c - 2to - \xr 

J-T/2 



l-T/2 

X e (ri - to) e (t2 - to) 
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Appendix A: Relative normalization 

Normalization requires evaluation of 
1 V 

v '''' ^ ' 



dkk^ 







(uj - eY + r2 



Since Wc ~ 2|kr| according to ([5]), we want u = 2k. 
Making a change of variables to rj = 2k — E with limits 
—oo to oo and selecting a contour that encloses the pole 
at 77 = —iT with T <^ E gives 



If T > the limits can be extended to ±00 and the 
0-functions imply that 



/3 = exp [-F (2tc - |x^|)] f dto exp [2rto] 

J —00 

= {2Ty' exp [-r {2t, - 2t< - |x,|)] 



where r> (t< ) is the larger (smaller) of ti and T2 . Since 
according to (H]) and 



2tc - 2t< = t> + t< - 2<< + 2 |x< I 
= f>-t<+2|x<|, 



II 



V 

2^ 



27ri 
4if 



VE^ 
167rr' 



Appendix B: Relative k- to x-space integral 
To evaluate (IT6l) we need 



13 = (2r) ^ exp [-r (i> - i< - |x^ I -f 2 |x< | 



12 



V^E^ 



d-^k 



exp (zk • r — iujt) 



IGnT 27r 
E^ ir 
X exp (—i2kt) . 

-- (exp 
47r r 1 



dkk 



uj-E + iV 
exp {ikr) — exp {^ikr) 



2k 



iT 



{tE + T/2){t--r 



— exp 



{iE + r)[t + -r 



Eq. ([23|) gives \xr \ = |x> | -I- |x< | , that is the distance be- 
tween the detectors equals the sum of the source-detector 
distances, so that 



/3 = (2r) ^ exp [-F {ty - t< - x> + |x< | 
= (2r)-iexp[-F(r>-r<)]. 
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